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We derived equations of motion corresponding to Bianchi-I cosmological models in the Lovelock 
gravity. Equations derived in the general case, without any specific ansatz for any number of spatial 
dimensions and any order of the Lovelock correction. We also analyzed the equations of motion solely 
taking into account the highest-order correction and described the drastic difference between the 
cases with odd and even numbers of spatial dimensions. For power-law ansatz we derived conditions 
for Kasner and generalized Milne regimes for the model considered. Finally, we discuss the possible 
influence of matter in the form of perfect fluid on the solutions obtained. 
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q . I. INTRODUCTION 

<n ; 

The idea that our space-time has more than four dimensions is not new in science. And since its first introduction 
by Kaluza [l[ and Klein Q it went a long way with different inspirations. Nowadays, the interest to such theories 
is mostly motivated by the superstrings and supergravity theories (see, e.g., Q). The presence of extra dimensions 
0^ ' makes it possible for superstrings theory to exist, but it leaves us with a gap between the number of dimensions the 
superstrings exist in, and our observable four-dimensional space-time. One of the possible ways to solve this problem 

■ is to make additional dimensions very small on cosmological scales by means of the cosmological theory. 

It is well-known that the Einstein tensor is, in any dimension, the only symmetric and conserved tensor 

O h . depending only on the metric and its first and second derivatives (with a linear dependence on second derivatives). 

If one drops the condition of linear dependence on second derivatives, one can obtain the most general tensor which 
bJ). satisfies other mentioned conditions - Lovelock tensor [7]. The corresponding Lagrangian (Lovelock Lagrangian), 
which generalizes the Einstein-Hilbert one, can be obtained by the same way - the Euler variation with respect to 
^v^j , the metric. The Lovelock tensor differs from the Einstein tensor only if space-time has more than four dimensions, 

■ so the Lovelock gravity can be considered as a natural generalization of General Relativity in a higher-dimensional 
t-H ' space-time. It contains quadratic and higher-order corrections in the Lagrangian, depending on the number of spatial 

■ dimensions. 

In the cosmological context, the Lovelock gravity has as well been intensively studied (see, e.g., [8-12]). In our 
paper we study the Bianchi-I-type cosmological model in Lovelock gravity. As we pointed out above, the anisotropic 
cosmological models can possibly solve the issue with small extra dimensions. But to construct a self-consistent model, 
freed from singularities and other obstacles, one needs to fully understand the dynamics of such a system. 

We decided to investigate Bianchi-I cosmological models. These are models with a simplest geometry and they can 
be used to understand common features of the model considered. The case with most general kind of geometry is 
t> • more interesting, of course, but is much more complicated as well. And knowing some features from the Bianchi-I 
model analysis could be a great asset while studying more general cases. Not to mention, exploring Bianchi-I models 
is profitable by itself - indeed, according to nowadays cosmological point of view, our Universe is spatially flat, so as 
Bianchi-I models are, while general case, which usually involves products of Lie algebras, is not. 

As a part of the analysis we demonstrate the difference between the cases with odd and even numbers of spatial 
dimensions. In our previous paper dedicated to the numerical investigation of the Gauss-Bonnet cosmologies |l3[ we 
noticed that there is a clear difference in the dynamics of (4+1)- and (5+l)-dimensional models. Now we demonstrate 
that this difference is generic for all orders of Lovelock correction and for any number of dimensions. 

Finally we take the matter in form of the perfect fluid into analysis. This decision is also motivated by our previous 
studies. In [l4[ we discovered the existence of "special solution" for pressureless fluid in (4+l)-dimensional Gauss- 
Bonnet cosmology; in this paper, among other things, we demonstrate that this is the common feature of models with 
even number of spatial dimensions. 

The structure of a manuscript is as follows: first, we derive the Lagrangian for the model with only one highest 
possible n taken into account as well as equations of motion. Then, we consider cases with even and odd numbers of 
spatial dimensions separately, demonstrate the difference in the structure of the equations and, as a consequence, a 
difference in their dynamics. In section IV we apply so-called power-law anzatz to our equations and obtain conditions 
for Kasner and generalized Milne regimes; conditions obtained generalize known (for low n) ones. Then we take into 
account matter in form of a perfect fluid and investigate the influence of it onto the dynamics of the system. In 
Discussion we discuss the topics relevant to the results obtained and finally in the Conclusion sum up all the results. 
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II. EQUATIONS OF MOTION 



Lovelock invariants have the form 



_fJlJ2...l2n W 1 ! 2 

2" hh - hn ilil 



..m 



(i) 



where Sj\ l j 2 '"J™ is the generalized Kronecker delta of the order 2n. The Lagrangian density have a form 



C 



-9 



^ a n L n , 



(2) 



where g is the determinant of metric tensor, a n is a constant of the order of Planck length in n dimensions and 
summation over all n in consideration is assumed. In our analysis we take into account only one highest possible 
correction - we are interested in the behavior in the vicinity of the cosmological singularity. In that regime higher 
curvature corrections are dominating over lower ones so one can consider highest possible correction only. In that 
case ([2]) contains only one term with highest possible n: 



C 



-gL r , 



(3) 



since we have only one n, we can renormalize L n to set a n to unity. As we mentioned in Introduction, we work with 
Bianchi-I-type metric in (D + l)-dimensional space-time: 



ds 2 



diag(-l, ai(t), a 2 2 {t),...a 2 D (t)), 



where di(t) is a scale factor, corresponding to i-th spatial dimension. From this metric one can build the following 
non-zero Riemann tensor components: 



Rom = -ai{t)'di(t), R t jij = ai(t)ai(t)aj(t)aj(t), 



(4) 



where Latin indices correspond to the spatial coordinates, the dot above corresponds to a derivation with respect to 
time and standard relations between indices for Riemann tensor are kept in mind (see, e.g., [HI). By raising indices 
in Eq. Q one can get 



0i ~ o,(t)' 



R ij = ai(t)aj(t) 



(5) 



Since we are working in a flat space-time, it is convenient to rewrite Eq. ([S]) in terms of the Hubble parameter 
Hi = ai(t)/ai(t): 



jdOi TJ I tt2 

Mm — Ui+Uii 



R lj - HiHj. 



(6) 



Now we can substitute into ([TJ , imply the properties of the generalized Kronecker delta and obtain the expression 
for n-th order Lovelock invariant for Bianchi-I-type metrics: 



L 



D 



D 



(2n-3)!!£ 



(Hi + H?) £ 

31>j2>...j2n-a 



. . Hi. 



(7) 



+(2n-l)!! J2 H kl H k2 ...H k2 

k 1 >k2>...k 2 r l 



where (2n — 1)!! is bifactorial of (2n — 1): (2n — 1)!! = (2n — 1) x (2n — 3) x • ■ • . There D is a number of spatial 
dimensions and D ^ 2n. Since we deal with only one (highest possible) n, we can neglect the l/2 n multiplier as well 
as cancel (2n — 3)!! with a reduction of (2n — 1)!! and Eq. (J7J become 
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-D 



D 

= E {h 

i=l 

+(2n-l) 



+ H?) E 

{j'l,j2: ■•■ j2n-2}^« 
Jl>j2>---32n-2 

E ^fel ^2 ■ • ' ^2 

kl >k 2 >...k2n 



H 3i H h 



. H, 



(8) 



The first term in Eqs. (J7]) and ([8]) comes from the products of a R^R 1 



products like R^Rfi 



■jk 

Since the same index can not appear twice, the products of a R^\R^ 



kind and the second one - from the 
U j ■ ■ . kind are impossible. 



Then one can substitute Eq. (|8|) into Eq. ((3]) and infer the equations of motion via varying; thus, the constraint 
equation takes a form 



(2n-l) H kl H k2 ...H k2n =0; 

k±>k2> ...k2n 

the dynamical equation corresponding to i-th spatial coordinate becomes 



(9) 



D 

E 

m— 1 



{Hm 
-(2n 



~ H m) E 

{j'l,j2, ••• j2n-2}^{i,m} 
jl>h>—hn-2 
1) E Hk 1 Hk 2 ■ ■ ■ H k . 

{k 1 ,k 2 ,...k 2 ,i}^i 
k 1 >k2>...k2 rl 



Hi. 



+ 



(10) 



III. THE DIFFERENCE BETWEEN THE D = 2n AND D = 2n + 1 CASES 



Eqs. ([9]) and (fT0|) might look complicated but in fact they are quite simple; to illustrate it let us rewrite them for odd 
and even numbers of spatial dimensions. For this one can notice, that any given n gives a maximal-order correction 
for an odd D = 2n + 1 and even D = 2n number of spatial dimensions. So for odd (D — 2n + 1) number of spatial 
dimensions, the constraint equation takes a form 



Y, H kl H k2 ...H k2n =0: (11) 

k 1 >k2>...k2n 

this is just the sum of all possible multiplications put in order with 2n multipliers constructed from Hubble parameters. 
One can see there are only (2n + 1) terms in this sum, and we omitted the (2n + 1) multiplier in front of sum in © 
since it is always non-zero. The dynamical equation corresponding to i-th spatial coordinate become 



D 

E 

m—1 



{jl,32, ■•• J2n-2}#{«,m} 
Jl>j2>-..j2n-2 

-(2n-l) H kl H k2 ...H k2n =0. 



Hi. 



(12) 



For the case considered, the second term in Eq. (fTU|) has only one possible combination - indeed, we have (2n+l) spatial 
dimensions, one of them is "forbidden" (since this equation corresponds to the i-th coordinate) and the remaining 2n 
coordinates have only one multiplication combination with 2n multipliers. 

For an even number of spatial dimensions (D = 2n) we have a constraint equation in the following form: 



H kl H k2 . . . H k . 2n — 0. (13) 

The same with the Eq. (jlip reasoning, we have only one multiplication with 2n multipliers constructed from Hubble 
parameters for this case. And this multiplication is equal to zero. This leaves us the first drastic difference between 
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the cases with odd and even numbers of spatial dimensions - the latter should have one of Hubble parameters always 
equal to zero. The dynamical equations for the odd-number case: 

D 

(H m + Hi) H h H j2 . . . H hn _ 2 = 0. (14) 

m— 1 

One can clearly see the difference between Eqs. (fTU)) and (TT2"1) : the last term vanished due to the absence of multipli- 
cations of this type and the r.h.s. of the first term is a multiplication with 2n — 2 terms: one index is "forbidden" 
because it's an i-th dynamical equation, another one comes into (H m + H^) part, consequently we are left exactly 
with 2n — 2 "free" Hubble parameters which form only one put in order multiplication with 2n — 2 multipliers. 

Therefore one can clearly see the difference between the cases with odd and even number of spatial dimensions - 
they differ both in the constraint equation as well as in dynamical ones. In [l4[ (and previously in [l3j]) we showed 
the difference between the (4+ 1) and (5 + 1) Gauss-Bonnet-dominated (L2) Universe; now we see that this difference 
is global and that it is applied to all orders of Lovelock gravity. 



IV. POWER-LAW ANZATZ 



Since the Kasner's paper [16j . the power-law anzatz for Bianchi-I model is widely used for various reasons. It implies 
the scale factors to be written as ai(t) — t Pi , so the Hubble parameters would be Hi — pi/t and the derivation of 
Hubble parameters with respect to time would have the form Hi = —pi/t 2 . Let us rewrite the equations of motion 
using this anzatz and obtain some important results. Thus, constraint equation of our system takes a form 



k 1 >k 2 >...k2n 



■Pk 2 



0; 



(15) 



we omit the (2n — l)/t 2n multiplier since only one n (the highest possible) is taken into account and we deal with a 
vacuum case. The dynamical equations become (again, we write down only one z-th equation) 



D 

E 

m— 1 



{p 2 m~Pm) E PnPh---Pn 

{jl,33, ■■■ J2r 1 -2}#{i,m} 

31>32>---hn-2 

+(2n-l) E Pk t Pk 2 ■ ■ • Pk 2n = 0; 

{k 1 ,k 2 .,...k 2n }^i 
ki>k 2 >...k 2n 



(16) 



and we omit the t 2pm ~ 2n multiplier. Needless to say, all the discussion in the previous section about the difference 
between the cases with odd and even numbers of spatial dimensions is valid for this anzatz, and the analogs of 
Eqs. (TTTI FT¥)) could be easily written. 

Now let us find the conditions for the analogs of Kasner regime. In finding them we followed [T3]- Let us consider 
the sum of all dynamical equations - the sum of all Eqs. (|16[) for all i. We will have a sum 



D 

E E 

[i] m=l 



(Pm Pro) E P31P32 ■ ■ ■ Pj 2n - 2 

{jl,32, ■■■ 32 n - 2 }^{i-m} 
3l>h>-hn~2 

(2n-l)E E Pk t Pk 2 ■ ■ -Pk 2n = 0. 

[i] {/Sl,/E2,...fe2n}#i 

ki>k2>...k2n 



(17) 



Taking (| 15[) into account, one can see that the second term in (jTTJ) is zero. As for the first term, it can be decomposed 
into two parts - with only linear multipliers and with only quadratic ones: 



(D + l-2n)J2 Pi E PhPh ■ ■ ■ Phn-2- 

i=1 {jl, 32, ■■■ J2n-2}^i 
il>i2>-i2r,-2 

-(2n- 1)(D + 1 - 2n) £ PiPhPj 2 ■ ■ ■ Pj 2n -2 = °; 

i>3l>32>---hn.-2 
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the coefficients before the terms are the result of simple combinatorics. Cancelling {D + 1 — 2n) one would obtain a 
relation, which depends only on the order of Lovelock correction and not on the number of spatial dimensions: 



D 



Y PllPh ■■■Pj2n-2 ~ (2/1- 1) Y PhPh--Phn-l =0 - ( 18 ) 



i=1 {jl, 32, ■■■ J2n-2}/i jl>32> ■■■32n-l 

jl>j2>—hn-2 

Now let's consider a sum 

D 

Y Pi Y PhPh-'-Ph n -i- 

1=1 J'l>j2>---J2r(-1 

This sum falls into the terms of two kinds - if i = {ji, J2, • ■ • J2n-i}, its terms are like the first term in (|18[) . if 
* 7^ {iii J2, • ■ • J2n-i} ~~ like the second ones. After some simple combinatorics, one can rewrite this sum as 

D D 

E Pi E PhPh ■ ■ -Phn-i =EP« E PjiPia • • •Pj 3 „-3 + 

*=1 il>j2>-..j2n-l i=l {jl, j 2 , ■■■ j2n-2>#« 

jl>32>— hn-2 

+2n PiPjx---Phn- 

i>'3l>32>---32n 

From (|15p it is obvious that the second term is equal to zero. So finally 

D D 

Y Pi Y PhPh---Phn-i=Y, p2 i Y PhPh---Phn-2- (19) 

»'=1 Jl>j2>..j2n- 

Combining (fT5|) and (TH?)) results in 



»'=1 3l>32>---32n-l J=l {il,j2, ■■■ 32»-2>#i 

J1>J2>— j2n-2 



■ ■ ■ p ^„-i = ( 2n - !) 51 ■ ■ ■ ' 

«=1 3l>j2>...j2n-l Jl>j2>---J2,,-1 



so cither 



^ Pi = (2n-1), (20) 



i=l 



or 



Y PhPh---Phn-i=°- ( 21 ) 



Jl>j2>...j2n-1 



First of the two is the condition for the analog of Kasner regime in (D + l)-dimensional Universe with Lovelock 
gravity of the order n, while the second corresponds to the generalized Milne solution (see the Discussion for more 
details). 



V. INFLUENCE OF MATTER 



In this section we will investigate the influence the matter in form of a perfect fluid exerts upon the dynamics of 
the system. The equations of motion in this case take the form 



6 



Guv 



8nG, 



(22) 



where Gfl is a D-dimensional gravitation constant, T^ v is the stress-energy tensor of a perfect fluid 

Tfj.u = (p + p)u^u v — pg^ v with an equation of state p — wp, (23) 

and G^v is a generalization of the Einstein tensor in Lovelock gravity (Lovelock tensor; for issues about its properties 
as a generalization of Einstein tensor, see, e.g., [Hj]): 



8g^ v 2 

So, the equations of motion would take a form: 
for constraint 



SH 1 

G^y = —g^yC where L is the Lagrangian ([3]). 



(24) 



(2n-l) ^ H kl H k2 . . . H k2n = p; 

ki>k2> ...fe2n 

and for dynamical equations (again, i-th equation) 



(25) 



E 

m— 1 



(H, 
+(2n-l) 



#m) E H 3i H h ■ ■ ■ H hn- 

{h,h, ■•• J2n-2}#{«,m} 
3l>32>...J2n.-2 

E H kl H k2 ...H k2n +p = 0. 

{k 1 ,k 2 ,...k 2n }^i 
ki>k2>...k 2 n 



(26) 



In the equations (1251 126[) . we renormalized the density p in order to avoid the 8ttGd/3 term - we are interested in 
the qualitative description so wc can do this. And of course, in addition to (|25l |26|) one needs a continuity equation: 



P+(P + P)E^ = - 



(27) 



The system (I25H27I) . along with initial conditions, completely determines the dynamics of the system and its past 
and future evolution. 

Let us rewrite (l25l [26| for the power-law anzatz. It's easy to see that Lovelock tensor scales as 



Goo = (2n-l)t- 2n J2 lPl 

%=* 2p '- 2n E 2 Pi, 



(28) 



where EiP» an< i E2P* are ^ ne corresponding sums in l.h.s. of Eqs. (|15p and (|16l) . respectively (one can easily verify 
that both terms in Eq. (|16|) scale the same). Solving (|2"T1) in the chosen anzatz, one can write down the scaling of the 
stress-energy tensor components: 



Matching (gSJ) and 



Tjj = wp t 2 Pi-( 1+w ^P\ 
one can find the condition under which system can be solved: 

Pi 



1 + w 



(29) 



(30) 



Finally, matching (|20p with (|30l) . one can find the "critical" value for the equation of state, with separate Lovelock 
and matter-dominating regimes (see the Discussion for more details): 



1 

2n - 1 ' 



(31) 
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VI. DISCUSSION 



We have derived the vacuum equations of motion I10[) for Bianchi-I-type metric for any number of spatial 
dimensions (D ^ 4) in the Lovelock gravity of any order n ^ 2. The corresponding Lagrangian for L 2 was first 
derived by Lanczos [l9j and for L 3 in Q. One can verify, that for the chosen metric our result coincides with the 
previously obtained ones. 

If one considers only the leading (in \ jt order) terms, then one needs to take into consideration only the highest- 
order Lovelock contribution, it would be n = [D/2]. Therefore we rewrote our equations of motion for the cases of 
odd D = In + 1 ri2"j) and even D = 2n (fT3l [T4"f numbers of spatial dimensions. One can easily see the difference 
between the two: while there are no restrictions from constraint on Hubble functions in case of odd numbers [of 
spatial dimensions] (Eq. (fTTj) ). the even- numbers case (Eq. (|13|)) sets a restriction: one of Hubble parameters should 
be always equal to zero. The dynamical equations also differ in structure ( (fT2|) vs. ([T4|) ). Our numerical investigation 
of (4+1) and (5+1) cases in Li [14| proved that there is a significant difference between two cases. We believe the 
difference in Li is due to the mentioned above reasons, so one should expect similar features from higher-dimension 
cases as well. 

When it comes to the anisotropic cosmological models in Lovelock gravity, two main cases are considered: product 
spaces and Bianchi-I models. As we claimed in Introduction, we consider the latter one. And for Bianchi-I-type 
models most of results are obtained in power-law anzatz. We also applied this anzatz to our solution and investigated 
it a bit. It has a form (|15[ I16p and coincides with previously found results for corresponding n (in addition to the 
previously mentioned papers, let us add [lfl| with a result for [D — 9, L4}). 

For the power-law anzatz we also found an expansion rate in terms of Kasner exponents: ^2 i p% = (2n — 1). This 
result also coincides with those previously reported for low n. Finally, we found the condition on pi, which corresponds 
to the generalized Milne solution. It is determined from two conditions - ([TBI and (|2"Tj) . And here we can again see 
the difference between the cases with odd and even numbers of spatial dimensions. Indeed, for D = 2n, the Eq. (fT5)) 
reduces to only one term directly implying that one of pi should be zero. With this, Eq. (f2~Tj) also reduces to only one 
term implying another pi should also be zero; the rest of pi are arbitrary. Thus, the resulting pi layout is as follows: 
(ax, ai, . . . , od_2, 0, 0). For D = 2n + 1, the situation is a bit different - since D = 2n + 1, we have D terms in 
Eq. (|15p . and nullifying one of p, will leave us with only one term in l.h.s. of Eq. (fT5"|) . and this eventually nullifies 
another pi. But now with two of pi set to zero, the l.h.s. of Eq. (|21[) is also left with only one term and this sets third 
of pi to zero. Therefore for odd numbers of spatial dimensions the resulting pi layout is as follows: (ax, a-i, . ■■ , 0,0-3, 
0, 0, 0). Finally, in a D- independent form: (01, ai, . . . , <Z2n-2) 0, . . . ). A proof that there are no solutions except the 
two mentioned above (Kasner-like with (|2T)|) and generalized Milne with (|2"Tjl ). is quite similar to the proof for L\ 13]. 

Finally we studied possible influence of matter in the form of perfect fluid on the dynamics of the system. We wrote 
the equations of motion ([2"5H2"?| and deepened the difference between the cases with odd and even numbers of spatial 
dimensions. In the presence of matter, even-dimensional cases change drastically - without the matter one of Hubble 
functions should be zero (see Eq. (|13p ). while even a small amount of matter changes the situation - in the presence 
of matter (p ^ 0) none of Hubble parameters can be zero (see Eq. (|25p). 

We as well found a condition for the system to be solvable in the power-law anzatz (Eq. (|3T))) ). Matching (f2T)|) and 
(13"0|. o ne can find a "critical" value for the equation of state, with separates Lovelock and matter-dominating regimes. 
In [lj] we demonstrated that this value separates two different asymptotic [leading to initial singularity] regimes - 
if w > w cr the initial singularity is isotropic (and the matter terms are dominating), and if w < w cr the initial 
singularity is anisotropic (and the Gauss-Bonnet terms are dominating). We believe that the same situation holds in 
the higher orders of Lovelock gravity, though some numerical simulations might be an asset. Finally, if one rewrite 
Eqs. (1231) and (|26l) for D = 2n and w = in terms of pi, they would become PkiPk 2 ■ ■ - Pk 2n = P f° r the constraint and 



for the dynamic equations (there we omit the powers of t). One can see, that the "special solution" pi = 1, which was 
described in [l4| for L\, is valid for the general case of Lovelock gravity in even numbers of spatial dimensions. 

And the last, but nonetheless a very important note. All through the present paper we claim that our results 
are obtained with only one highest-order Lovelock correction taken into account. In reality, some of them can be 
generalized for a mixture of Lovelock corrections with different n. If one wants to deal with some different Lovelock 
corrections at once (this would be the case if one wants to study the dynamics not only in the vicinity of an initial 
singularity), one should use ^ instead of © as a Lagrangian density and substitute all individual L„ needed into 
it. But in this case one need to use (J7]) expression for L n since we cancel 2~™ x (2n — 3)!! terms in © in order to 



D 




m—1 



8 



simplify it for the case where only one n is used. Now if we combine this n-dependent multiplier with a n into new 
constant a n , equations take a form 



where Pi are the equations of motion corresponding to Li - l.h.s. of Eq. © for a constraint or Eq. (ITU)) for a dynamical 
equation. Thus, equations can be easily constructed from the already derived ones. The same comes for the matter- 
filled universe, but in this case in the r.h.s. of ([52)1 would be p for a constraint or (— p) for the dynamical equations 
(compare with $?5§ and (f2"6")0 . 

From the way we derived the equations for the power-law anzatz one can see that the power-law analogs of equations 
cannot be generalized for the mixture of n. This comes out from the way we derived them - as you remember, we 
omitted the powers of t since they were the only ones for vacuum case or the same with those in r.h.s. for the matter 
case. But with a mixture of n we would have different powers of t for different n and so we cannot cancel them 
anymore. With this, the results obtained from the power-law anzatz (like different sums of pi), cannot be applied for 
the case with mixture of n. Yet, the results for the highest possible n can be used as an asymptotic regime for this 
case. 



In this paper we investigated Bianchi-I cosmological models in Lovelock gravity. We were first to derive the 
Lagrangian and the equations of motion for the model considered: Bianchi-I model in Lovelock gravity with only 
highest possible correction taken into account. This by itself is a great achievement - previously, when one deals with 
Bianchi-I models in Lovelock gravity, only one specific order was involved. This is linked with a difficulty in obtaining 
equations from the most general case. In our paper we imply symmetries of Bianchi-I model to obtain our results; 
our result is correct only for Bianci-I case but it cover all possible n as well as is valid for any number of spatial 
dimensions. 

We analyzed obtained equations of motion separately for the cases of odd and even number of spatial dimensions 
and proved that there is a difference between these two cases. 

After applying power-law anzatz to our equations of motion we were able to obtain conditions for the Kasner as 
well as generalized Milne regimes. Results obtained matches those for low-n cases and generalize them. 

Finally, we took the matter in form of perfect fluid into consideration. For this case we also obtained equations 
of motion; condition this equations can be resolved under and the "critical" value for the equation of state which 
separates Lovelock-dominating from matter-dominating regimes. Results of the analysis of the model with matter 
make the mentioned above difference between D = 2n and D = 2n + 1 cases even deeper. 

We also discussed the way of constructing the equations of motion for a mixture of n from the obtained equations 
both for vacuum and for matter cases. 

All results obtained matches with existing results for low-n cases and generalize them for all possible n. 



I would like to thank A.V. Toporensky and I.V. Kirnos for a stimulating discussion and S.V. Karpov for technical 
help. 
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